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Abstract. 

Within the nuclear density functional theory (DFT) we study the effect of reflection- 
asymmetric shapes on ground-state binding energies and binding energy differences. To this 
end, we developed the new DFT solver AXIALHFB that uses an approximate second-order 
gradient to solve the Hartree-Fock-Bogoliubov equations of superconducting DFT with the 
quasi-local Skyrme energy density functionals. Illustrative calculations are carried out for even- 
even isotopes of radium and thorium. 



1. Introduction 

The goal of low-energy nuclear theory is to provide a comprehensive theoretical framework to 
explain properties of atomic nuclei and their reactions at the nucleonic level. The roadmap 
towards this goal involves three theoretical strategies: ab initio methods, configuration- 
interaction techniques, and the nuclear DFT ^ E|- Due to exploding dimensions of the 
configuration space, properties of complex heavy nuclei with many valence particles are best 
described by the nuclear DFT and its various extensions 0. The main idea of DFT is to 
describe an interacting system of fermions via its densities rather than the many-body wave 
function. The energy of the many body system can be written as a density functional, and 
the ground state energy is obtained through the variational procedure. Figure [T] illustrates the 
strategy behind the nuclear DFT. A major thrust in this area is to enhance the predictive power 
of DFT by firmly connecting the nuclear energy functional (EDF) to the nuclear interactions and 
to optimize the low-energy coupling constants of the microscopically based EDF by comparing 
results of theoretical simulations with a selected set of observables [H [5] . 

The nuclear DFT is based on the self-consistent mean-field approach rooted in the self- 
consistent Hartree-Fock (HF) or Hartree-Fock-Bogoliubov (HFB) problem. The self-consistent 
HFB equations allow us to compute the nuclear ground state and a set of elementary quasi- 
particle excitations that can be used to build excitations of the system. The HFB equations 
constitute a system of coupled integro-differential equations that can be written in a matrix form 
as an eigenvalue problem, where the dependence of the HFB matrix elements on the eigenvectors 
induces nonlinearities. 
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Figure 1. The nuclear DFT strategy diagram. 



Since all the nucleons are treated on equal footing and contribute to single-particle and 
pairing mean fields, a few hundred nucleonic wave functions have to be considered which makes 
systematic calculations quite demanding. We have at our disposal a suite of optimized HFB 
solvers working in various geometries [6j. All our DFT solvers adopt an iterative approach 
to solving the HFB equations [3 [HI O [TOl [Til [12 • Nuclear configurations (characterized by 
shape deformations, angular momentum, temperature, single-particle occupations, etc.) can be 
selected by means of external constraints [13] • By employing modern DFT infrastructure global 
nuclear properties can be rapidly computed across the nuclear landscape [14j. 

One of the fundamental properties of an atomic nucleus is its shape. The DFT description 
of nuclei is performed in the frame of reference of the nucleus, the intrinsic frame, in which the 
nucleus may acquire a deformed shape. The concept of shape deformation is ultimately related to 
the spontaneous symmetry breaking effect known in many areas of physics. Symmetry breaking 
solutions may appear variationally when, in a mean-field configuration respecting the original 
symmetries of the nuclear Hamiltonian, degenerate single-particle orbits are strongly coupled to 
collective vibrations ^,[1611171118]. In heavy nuclei, the most important collective modes are the 
quadrupole and octupole vibrations, leading to intrinsic mean fields with non-zero quadrupole 
and octupole moments. If a nucleus has nonzero quadrupole moments, its non-spherical shape 
can be characterized by ellipsoidal deformations [20]. The presence of octupole moments 
is indicative of reflection-asymmetric deformations [2T]. While quadrupole deformations are 
common in ground states across the mass table, octupole deformations are more concentrated 
in particular regions of the chart of the nuclides [211 [22] • 

This paper introduces a new DFT solver AXIALHFB which solves the HFB problem by means 
of an approximate second-order gradient method [23j. Using AXIALHFB with Skyrme EDFs, we 
study the effect of reflection-asymmetric moments on nuclear binding energy in even-even Ra 



and Th nuclei. These nuclei were chosen for their well known strong octupole correlations. The 
binding energies obtained in the nuclear DFT are further used to compute mass filters such as 
two-nucleon separation energies and the double-difference indicator SVpn [Ml ESI '26] . 

This paper is structured as follows. In Sec. [2] we give a brief overview of HFB equations, 
introduce the new DFT solver AXIALHFB, and provide necessary details about our DFT 
calculations. The results for the mass filters are presented in Sec. [3l Finally, conclusions and 
prospects for future work are given in Sec. [H 

2. DFT Calculations 

The following is a succinct description of the HFB approach. For a more in-depth discussion, the 
reader is referred to Refs. [3l|20j. The goal is to find the ground state energy of a nucleus (A^, Z) 
subject to constraints. The energy functional is a three-dimensional spatial integral of local 
energy density T-L that is a real, scalar, time-even, and isoscalar function of local densities and 
their derivatives. The energy density consists of a kinetic term, the Skyrme energy functional 
representing the effective nuclear interaction between nucleons, and the Coulomb term. A 
Skyrme functional depends on a number of local densities: nucleonic densities, kinetic densities, 
spin densities, spin-kinetic densities, current densities, tensor-kinetic densities, and spin-current 
densities [27j. Since pairing correlations are considered, the EDF also includes pairing terms 
defined in terms of pairing densities. Variation of the binding energy leads to the HFB equations: 



where h is the HF Hamiltonian containing the self-consistent field, A is the pairing field, A is the 
Fermi level, Ej^ are one-quasiparticle energy eigenstates, and (Uk^Vk) are two-component HFB 
eigenvectors. To restore the particle symmetry spontaneously broken in the HFB approach we 
use the Lipkin-Nogami scheme [28l |29j implemented as in [30j . 

The present calculations were done using a parity-breaking axial DFT solver AXIALHFB 
recently developed by Stoitsov using the gradient method routine for solving the HFB equations 
^ provided by L.M. Robledo [23j. Like in the parity-conserving axial HFB code HFBTHO 
[8j, AXIALHFB expresses the HFB equations in a large configuration space of single-particle 
eigenstates of the deformed harmonic oscillator basis. Unlike HFBTHO, however, which employs 
the direct diagonalization technique, AXIALHFB uses an approximate second-order gradient 
method [23i i3i]- This method is particularly well suited to deal with multiple constraints 
[20j : hence, it does not require special techniques such as the augmented Lagrangian method 
[13] implemented in HFBTHO. 

As an individual constrained HFB problem can be solved on a single processor, the 
computation of a potential energy surface (PES) in a space of collective coordinates is 
embarrassingly parallel (see Fig. [2j). This also makes it possible to produce deformed HFB 
mass tables - including reflection asymmetric shapes - on leadership-class supercomputers such 
as JAGUAR at Oak Ridge National Laboratory's Leadership Computing Facility. 

The goal was to find the ground state (g.s.) energy, i.e. the global minimum on a 2D 
potential energy surface. To guarantee that a global, rather than local, minimum is found, 
one has to compute and scan the PES in two coordinates: the axial mass quadrupole moment 
Q20 (characterizing the elongation of the nucleus expressed in terms of deformation /32) and 
the mass octupole moment Qso (characterizing the mirror asymmetry, or pear-like deformation 
/3s)- Figure El shows three energy surfaces for different isotopes of Radium calculated using the 
Skyrme EDF SLy4 [32j together with a density dependent mixed pairing functional [33j. All of 
the calculations here were done using a large oscillator basis of 21 oscillator shells (2024 basis 
states). We performed 525 constrained calculations for different combinations of Q20 and Qso 
(indicated by white dots on the plots in Fig. [2|). The range of Qso-values was identical for 
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Figure 2. The contour maps of the 2D potential energy surfaces of ^^^Ra (left), ^^^Ra 
(middle), and ^^^Ra (right) calculated by constraining both the quadrupole and octupole 
moments characterized by shape deformations /32 and /Ss, respectively. The mesh points at 
which calculations were performed are marked by white dots. The energy (in MeV) is shown 
relative to the ground-state minimum. The Skyrme functional SLy4 [32j and mixed pairing [ 33] 
were used. 
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Figure 3. Nucleon density distributions (in cylindrical coordinates r and z) calculated for ^^^Ra 
(/32 = 0.00, 13s = 0.00), 226Ka (/32 = 0.20, = 0.198), and ^^ORa (/32 = 0.24, /Ja = 0.0) showing 
spherical, quadrupole and octupole deformed shapes. 

each they went from to 6000 fm^ in steps of 250 fm^. Following the minimization on the 
resulting grid, unconstrained calculations have been carried out starting from the grid minimum: 
in this way the precise total g.s. binding energy was obtained and used to compute mass filters. 

Figures [2] and [3] illustrate three situations typical to all nuclei considered [35[35]- In general, 
for a given isotopic chain, a transition is expected with increasing neutron number from spherical 
shape to quadrupole well-deformed shapes, via reflection-asymmetric minima. The PES for 
^^^Ra is characteristic of a spherical system, which is fairly soft in the quadrupole-octupole 
direction. For ^^^Ra, our calculations yield a well-developed minimum with stable quadrupole 
and octupole deformations. As seen in Fig. [3l in its ground state ^^^Ra is predicted to have a 
pear-like shape. The nucleus ^^^Ra is strongly deformed but reflection-symmetric. Its shape is 
indicative of appreciable quadrupole and hexadecapole g.s. moments. 



3. Results 

The objective of this work was to determine whether reflection- asymmetric g.s. moments may 
have an affect on binding energy differences (mass filters) of even-even Ra and Th isotopes. 




Figure 4. Results of HFB calculations using the EDF SLy4 [32j and mixed pairing [33j for the 
even-even isotopes of Ra and Th. Top: g.s. deformation parameters ^2 and obtained in the 
full minimization allowing reflection-asymmetric shapes (Ps 7^ 0) compared with those obtained 
by assuming mirror symmetry (Ps — 0). Bottom: 6Vpn for reflection-symmetric (triangles) and 
asymmetric (circles) shapes compared with experimental values (x) of Ref. [36j. 

Figure [4] (top) shows the equilibrium deformations /32 and /3s predicted with the Skyrme 
functional SLy4. It is seen that the octupole deformation maximizes around A^=138 and then 
gradually decreases with the neutron number. The effect of the octupole deformation on the 
quadrupole moment appears to be very small. This is evident by comparing the g.s. /32-values 
obtained in the full minimization (/33 7^ 0) with those obtained by assuming reflection symmetric 
shapes (/33 = 0). The bottom part of Fig. [H shows the filter SVpn- 
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which can be viewed as an approximation to the mixed partial derivative [26j: 

In general, the agreement with experiment for 6Vpn is improved when octupole correlations are 
considered. The spike around N = 130 in the f3s ^ results is due to the rapid transition 
between spherical and deformed shapes; this effect - typical to mean-field calculations - is 
supposed to be washed out if beyond-mean- field effects are taken into account. 

Finally, Fig. [5] shows predicted and experimental two-neutron and two-proton separation 
energies: 

S2n = BiZ,N)-BiZ,N-2), S2p = B{Z,N)-BiZ-2,N). (4) 
For these mass filters, the effect of octupole correlations is very small. 



SVpn = 7 [B{Z, N) - B{Z, N-2)-B{Z- 2, N) + B{Z-2,N- 2)] , (2) 
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Figure 5. Similar to Fig. [H except for two-neutron (top) and two-proton (bottom) separation 



energies. 



4. Conclusions 

The large-scale simulations utilizing leadership-class supercomputers have transformed the 
DFT-based description of nuclei. A close coupling of math and computer scientists with 
physicists, focused on a close coupling with present and future experiments and observations, 
has dramatically advanced the field [37j. 

This work presents selected results of large-scale DFT calculations including reflection- 
asymmetric shape degrees of freedom. We focus on the influence of parity-breaking intrinsic 
moments on the binding energy differences. Our calculations suggest that SVpn is affected by 
octupole effects, bringing it closer to experiment, whereas S2n and S2p are not. Future work will 
involve performing similar calculations across the nuclear chart, also for odd- A nuclei. This will 
require extending AXIALHFB to the odd-particle case [38]. We also intend to study the model 
dependence of our predictions by considering different energy density functionals. 
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